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We report results for K/tt fluctuations from Au+Au collisions at y / ijvjv = 19.6, 62.4, 130, and 
200 GeV using the STAR detector at the Relativistic Heavy Ion Collider. Our results for K/n 
fluctuations in central collisions show little dependence on the incident energies studied and are on 
the same order as results observed by NA49 at the Super Proton Synchrotron in central Pb+Pb 
collisions at y / ijvjv = 12.3 and 17.3 GeV. We also report results for the collision centrality dependence 
of K/n fluctuations as well as results for K + /ir + , K~ /n~ , K + , and K~ /ir + fluctuations. We 
observe that the K/n fluctuations scale with the multiplicity density, dN/drj, rather than the number 
of participating nucleons. 
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PACS numbers: 25.75.-q, 24.60.Ky 



Strangeness enhancement has been predicted to be one 
of the important signatures of the formation of the quark 
gluon plasma (QGP) 0, & 0, i, H . Recently a maximum 
in the ratio of the yields of K + and 7r + , K + /tt + , has 
been observed in central Pb+Pb collisions near ^/sWn — 
7 GeV @. Dynamical fluctuations in the event-by-event 
K/tt ratio in central Pb+Pb collisions at energies near 
\J snn — 7 GeV are larger than those predicted by the 
transport model UrQMD using the observable <Jd yn 0- 
The K/tt ratio is defined as K/tt = (N K + +N K - )/(N 7T '+ + 
N n -), where (N K + + N K ) is the number of charged 
kaons in one event and (N^+ + N n -) is the number of 
charged pions in the same event. The observable <Jd yn Hi 
is defined as 

CTdyn = Sign (aj ata - CT^cd) ^kdata-^mixcdl i 1 ) 

where Odata is the relative width (standard deviation di- 
vided by the mean) of the K/tt distribution for the data 
and <7 m ixed is the relative width of the K/tt distribu- 
tion for mixed events. These observations have gener- 
ated speculation that a phase transition from hadronic 
matter to quark-gluon matter may be taking place at in- 
cident energies around ^snn = 7 GeV @. The study 
of dynamic fluctuations in the event-by-event K/tt ratio 
may produce information concerning QCD phase transi- 
tions such as the order of the transitions and the location 
of the transitions, and may lead to the observation of the 
critical point of QCD [I Gil. 

In this paper, we report results for dynamic fluctua- 
tions of the K/tt ratio in central Au+Au collisions at 
^/snn — 19-6, 62.4, 130, and 200 GeV using the quantity 
Cdyn @| ■ These results are compared with the results of 
NA49 for dynamical K/tt fluctuations in central Pb+Pb 
collisions To study the collision centrality dependence 
of K/tt fluctuations, we propose a new variable, Vdyn, Kir, 
which quantifies the deviation of the fluctuations in the 
number of pions and kaons from that expected from Pois- 
son statistics. This variable is defined as 

_ (N K (N K - 1)) (jV, (N w - 1)) 

dyn <^ ~ (N K f + (Nit) 2 ^ 

o (NicN n ) 
(N K ) (AM 

where N v is the number of pions in each event and Nk is 
the number of kaons in each event. The properties of Eq. 
2 are discussed at length in Ref. [ll|]. Negative values 
of Vdyn.KTT imply that the third term in Eq. 2 involving 
K — tt correlations dominates, while positive values of 
Vdyn,KTr imply that the first two terms involving the joint 
correlations K — K and tt — tt dominate. We present re- 
sults for the collision centrality dependence of Vd y n,Kn for 
Au+Au collisions at 62.4 and 200 GeV. To gain insight 
concerning the origins of these K/tt fluctuations [l2j , we 
also present the collision centrality dependence of fdyn, Ktt 
for K+/tt+, K-/tt-, K+/tt~, and K~/^ + . 



Au+Au collisions were studied at ^/snn — 19.6, 62.4, 
130, and 200 GeV using the STAR detector at the Rel- 
ativistic Heavy Ion Collider (RHIC). All data presented 
here were taken using a minimum bias trigger. Events 
accepted took place within ±15 cm of the center of the 
STAR detector in the beam direction. Collision central- 
ity was determined using the number of charged tracks 
within \r}\ < 0.5. Nine centrality bins were used corre- 
sponding to 0-5% (most central), 5-10%, 10-20%, 20-30%, 
30-40%, 40-50%, 50-60%, 60-70%, and 70-80% (most pe- 
ripheral) of the reaction cross section. To be able to plot 
our results versus dN/drj, we associate fully corrected val- 
ues for dN/drj from previously published work with each 
collision centrality bin [13j. For the 19.6 GeV and 130 
GeV data sets, only results from the most central bin are 
presented. All tracks were required to have originated 
within 3 cm of the measured event vertex. Only charged 
particle tracks having more than 15 space points along 
the trajectory were accepted. The ratio of reconstructed 
points to possible points along the track was required to 
be greater than 0.52 to avoid split tracks. Charged pions 
and charged kaons were identified using the specific en- 
ergy loss, dE/dx, along the track and the momentum, p, 
of the track. 

Charged pions and kaons were selected with transverse 
momentum 0.2 < p t < 0.6 GeV/c and pseudorapidity 
\rj\ < 1.0. Particle identification was accomplished by se- 
lecting particles whose specific energy losses were within 
two standard deviations of the energy loss predictions for 
a given particle type and momentum. Particle identifica- 
tion for pions (kaons) also included a condition that the 
specific energy loss should be more than two standard de- 
viations away from the loss predicted for a kaon (pion). 
In addition, electrons were excluded from the analysis 
for all cases. Particles were excluded as electrons if the 
specific energy losses were within one standard deviation 
of the energy loss predictions for electrons. We identify 
~90% of the pions in our acceptance. We identify ~50% 
of the kaons at pt = 0.2 GeV/c and ~75% of the kaons 
at p t = 0.6 GeV/c in our acceptance. We calculate that 
the fraction of pions resulting from misidentified kaons 
is negligible while the fraction of kaons resulting from 
misidentified pions is 6.5%. The electron cut did not af- 
fect the pions significantly, but excluded 25% of the kaons 
for the 200 GeV Au+Au case and 35% of the kaons for 
the 62.4 GeV Au+Au system. The remaining electron 
contamination is negligible. 

In Fig. [TJ we show the distribution of the event-by- 
event K/tt ratio for central Au+Au collisions (0-5%) at 
^snn = 200 GeV compared with the same quantity from 
mixed events. Mixed events were created by taking one 
track from different events to produce new events with 
the same multiplicity that have no correlations among 
particles in those events. Mixed events were produced 
using ten bins in collision centrality and five bins in event 
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FIG. 1: (Color online) The event-by-event K/n ratio for 
200 GeV Au+Au central collisions (0-5%) compared with the 
same quantity calculated from mixed events. The inset shows 
the ratio of the distribution from real events to that from 
mixed events. The errors shown are statistical. 
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FIG. 2: (Color online) Measured dynamical K/n fluctuations 
in terms of <7dyn for central collisions (0 - 5%) of 19.6, 62.4, 
130, and 200 GeV Au+Au compared with the central colli- 
sions (0 - 3.5%) of Pb+Pb from NA49 and the statistical 
hadronization (SH) model of Ref. The solid line rep- 

resents the relationship of the incident energy dependence of 
Cdyn in central collisions to the collision centrality dependence 
of fdyn,if7r at higher energies. Both statistical (vertical line 
with horizontal bar) and systematic (no vertical line) error 
bars are shown for the experimental data. 



vertex position. The distribution for the data is wider 
than the distribution for mixed events. Similar results 
were obtained at the other three incident energies; 19.6, 
62.4, and 130 GeV. The same analysis techniques are 
applied to the mixed events as are applied to the data. 

The results for <7dyn from central Au+Au collisions at 
19.6, 62.4, 130, and 200 GeV are shown in Fig. [U along 
with similar results for central Pb+Pb collisions at ^/snn 
= 6.3, 7.6, 8.8, 12.3, and 17.3 GeV from the NA49 Collab- 
oration [7[. Statistical and systematic errors are shown 
for both the NA49 results and the STAR results. The 
systematic errors for a^ yn are discussed in the presenta- 
tion of the results for ^d yn ,if7r- 
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FIG. 3: (Color online) Measured dynamical K/n fluctuations 
in terms of fdyn.A'ir for 62.4 and 200 GeV Au+Au compared 



with <Tdyn from central Pb+Pb collisions at 6.3, 7.6, 



12.3, 



and 17.3 GeV from NA49 0]. Statistical errors are shown for 
the STAR data. Statistical and systematic errors are shown 
for the NA49 results. The solid line corresponds to a fit to 
the STAR data of the form c + d/ (dN/dr)). 



In Fig. [51 we find that the NA49 results show a strong 
incident energy dependence while the STAR results show 
little dependence on the incident energy. The STAR re- 
sults are consistent with the highest energy NA49 re- 
sult, although the statistical error bar for the 19.6 GeV 
case is large. In this figure, we compare the statisti- 
cal hadronization model results of Torrieri [l4| to the 
experimental data. We see that when the light quark 
phase space occupancy, j 9 , is one, corresponding to equi- 
librium, the calculations underestimate the experimental 
results at all energies. When j q is varied to reproduce 
the excitation function of K + /n + yield ratios and the 
excitation function of temperature versus chemical equi- 
librium over the SPS and RHIC energy ranges 0, [l5[ , 
the statistical hadronization model correctly predicts the 
dynamical fluctuations at the higher energies but under- 
predicts the NA49 data at the lower energies, supporting 
the conclusion that the lower energy fluctuation data are 
anomalous [7]. In these fits, j q > 1 (chemically over- 
saturated) for ^/snn < 9 GeV and j q < 1 (chemically 
under-saturated) for ^snn > 9 GeV. 

We propose to study K/n fluctuations using a variable 
that does not involve the K/tt ratio directly. We choose 
to employ the variable fdyn,AT7r , which is similar to the 
observable V-\ ^dyn [3] used to study net charge fluctu- 
ations. Vdyn,KTT does not require mixed events and does 
not depend on detector efficiencies. In Fig. [3l we show 
Vdyn. Ktt for 62.4 and 200 GeV Au+Au collisions plotted 
as a function of dN/dr). We estimate the systematic error 
in Vdyn,K-K to be 15% due to losses from the electron cut. 
Using HIJING T7| , we estimate that the effect of feed 
down on t'dyn.A'Tr from weakly decaying particles is 9%. 
HIJING calculations show that increasing the accepted 
range in p t from 0.2 < p t < 0.6 GeV/c to 0.1 < p t < 2.0 
GeV/c decreases Vdyn^n by less than 5% at both 62.4 
and 200 GeV. 
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In Fig. [31 we plot the NA49 results using the identity 
CT dyn = v dyn.Kn- We verified the validity of this identity 
explicitly by calculating both quantities for 62.4 and 200 
GeV Au+Au collisions. We find that <J^ yn — Vdyn.K-rr 
within errors for all centrality bins except the two most 
peripheral bins at 62 GeV. In Fig. G2 we make the cor- 
respondence between the incident energy for the NA49 
results and dN/dr] using the systematics in Ref. [l8j . 
The solid line in Fig. [3] represents a fit to the STAR data 
of the form c + dj (dN/dr]) where c and d are constants. 
The fit for Vd yX L,K-K versus dN/dr] has a % 2 of 26.6 for 16 
degrees of freedom. If we make a similar fit for Vdyn,Kir 
versus A part , we obtain a \ 2 01 50.7 for 16 degrees of 
freedom. Thus, the fit for Vdyn,K-n versus dN/dr] is sig- 
nificantly better than the fit for Vdyn.K-K versus A par t- 
The NA49 results shown in Fig. [3] show a steeper depen- 
dence on dN/dr] than the STAR data and have a \ 2 of 
50.8 for five degrees of freedom compared to the best fit 
to the STAR data. 

Using the results of the fit for dN/dr] along with the 
systematics in Ref. |l8[ and the relationship cr^ = 
^dyn.ifTD we can draw the solid line shown in Fig. [3J The 
line relates the incident energy dependence of <7d yn in 
central collisions to the collision centrality dependence of 
ffjyn.ifir a t higher energies. We observe a slight increase 
in the systematic behavior of <7d yn based on scaled re- 
sults from STAR as the incident energy is lowered, while 
the NA49 results show a steeper increase as the energy 
is lowered. 

In order to gain insight into the origin of these K/n 
fluctuations, we calculate Vdyn,K-n for K + /n + , K~/n~, 
K + /n~, and K~ /n + . We observe that, within errors, 
^dyn,K+7r+ is equal to fdyn,K-7r- and 2/ dy n,K+7r- is equal 
to Vdyn, K-ir+ ■ We report the average of the fluctuations 
of the ratios K + /n + and K~ /n~ as same sign and the av- 
erage of the fluctuations of the ratios K + /n~ and K~ /n + 
as opposite sign in Fig. 2] along with results for K/n as 
a function of dN/dr]. Because Vdyn,K-n scales with the in- 
verse multiplicity, we multiply our results for t'dyn.ifTr by 
dN/dr] to study the collision centrality dependence more 
effectively. 

The scaled Vdyn,Kir results for all cases in Fig. 2] be- 
come more positive as the collisions become more cen- 
tral. The scaled results for the summed signs are al- 
ways positive and increase as the collisions become more 
central. The scaled results for the opposite sign are al- 
ways negative indicating a strong correlation for oppo- 
site sign particles. One might expect such negative op- 
posite sign correlations from processes such as the de- 
cay K* (892) — > K + + n~ . The scaled Vdyn,Kir for the 
same sign are slightly negative in peripheral collisions 
and slightly positive in central collisions, crossing zero 
around dN/dr] = 400. The fact that ^ dy n,-R"7r f° r same 
sign particles is close to zero indicates that the correla- 
tions between same sign particles are small. 

Also shown in Fig. [3] are filtered HIJING calculations 
for 62.4 and 200 GeV Au+Au collisions. These cal- 
culations include the acceptance cuts of \r)\ < 1.0 and 
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FIG. 4: (Color online) The dN/dr] scaled dynamical K/n fluc- 
tuations for summed charges (stars), same signs (circles), and 
opposite signs (squares) as a function of dN/dr]. The errors 
shown are statistical. The open and filled symbols refer to 
Au+Au collisions at 62.4 GeV and 200 GeV respectively. The 
dash-dot, dotted, and dashed lines represents HIJING calcu- 
lations for summed charges, same signs, and opposite signs 
respectively. 



0.2 < p t < 0.6 GeV/c. In contrast to the experimen- 
tal results, HIJING predicts little or no collision cen- 
trality dependence for (dN / d-q)vdyn,K-K ■ Because of this 
lack of centrality dependence, we show the average of 
the HIJING predictions as a horizontal dot-dashed line. 
HIJING predicts that the same sign fluctuations are al- 
ways zero (dotted line) and the opposite sign fluctuations 
are always negative (dashed line). HIJING does not in- 
clude resonances that decay to K + + n + or K~ + 7r~. 
One explanation for the increase in the measured same 
sign and opposite sign fluctuations scaled with dN/dr] 
may be that in peripheral collisions, products of the de- 
cay of resonances emerge without further interaction, 
leading to negative values of Vdyn,Kir ■ In central col- 
lisions, the daughters of the decay of resonances are 
rescattered in or out of our acceptance, leading to more 
positive values of Vdyn.Kn- For example, the decay of 
Ki (1270) + -> K+ + p° -> K+ + 7T+ + 7T" would lead to 
negative same sign fluctuations in peripheral collisions 
but not in central collisions. 

In conclusion, we find that the fluctuations in the K/n 
ratio for central Au+Au collisions at ^/s^n — 19.6, 62.4, 
130, and 200 GeV are of the same order as the fluctu- 
ations observed in central Pb+Pb collisions at ^/snn 
= 6.3, 7.6, 8.8, 12.3, and 17.3 GeV using the variable 
0dyn, but the Pb+Pb results show a stronger incident 
energy dependence. The statistical hadronization model 
of Ref. [lj] cannot reproduce the incident energy de- 
pendence of these fluctuations. The collision centrality 
dependence of K/n fluctuations for Au+Au collisions at 
\J sjviv = 62.4 and 200 GeV as characterized by the vari- 
able Udyn, Kit seems to scale with dN/dr]. Relating the 
observed centrality dependence of these fluctuations ob- 
served in Au+Au collisions at ^/sJTn — 62.4 and 200 



6 



GeV using Vdyn,Kn to the incident energy dependence 
of fluctuations in central collisions using <7d yn , we find 
a smooth scaling, decreasing slightly with increasing in- 
cident energy. The scaled values fall in the middle of 
the fluctuations observed by NA49 between ^Jsmn = 6.3 
and 17.3 GeV. dN/drj is a good measure of the freeze-out 
volume and we have shown that the fluctuations in K/tt 
scale with dN / drj. More measurements are required to 
demonstrate if there is any discontinuity in K/tt fluctu- 
ations as a function of incident energy. Vdyn.Kn results 
using pions and kaons with the same sign are close to 
zero while results using Vdyn,K-n for opposite sign pions 
and kaons are negative. The results for fdyn,K7r scaled 
by dN/ dr\ become more positive as the collisions become 
more central due to rescattering, unlike the predictions 
of H1JING that show no centrality dependence. These 
results may indicate that, due to later stage hadronic 
rescattering, the decay products of resonances are less 



likely to survive in central collisions than in peripheral 
collisions. 
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